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Superconductivity up to 30 K in charge neutrally doped BaFe2(Asi_a;Pa;)2 has been ascribed to 
chemical pressure, caused by the shrinking unit cell. But the latter induces no superconductivity in 
(Bai_a;Sr2,)Fe2As2 in spite of the same volume range. We show that the spin-density-wave (SDW) 
state of BaFe2As2 becomes suppressed in BaFe2(Asi_a;P2,)2 by a subtle reorganization of the crystal 
structure, where arsenic and phosphorus are located at different coordinates zab, zp. High-resolution 
X-ray diffraction experiments with BaFe2(Asi-a;Pi,)2 single crystals reveal almost unchanged Fe-P 
bonds, but a contraction of the Fe-As bonds, which remain nearly unchanged in (Bai_i:Sri)Fe2As2. 
Since the Fe-As bond length is a gauge for the magnetic moment, our results show why the SDW 
is suppressed by P-, but not by Sr-doping. Only the Fe-P interaction increases the width of the 
iron 3d bands, which destabilizes the magnetic SDW ground state. The simultaneous contraction 
of the Fe-As bonds is rather a consequence of the vanishing magnetism. Ordered structure models 
of BaFe2(Asi_i,Pa;)2 obtained by DFT calculations agree perfectly with the single-crystal X-ray 
structure determinations. The contraction of the Fe-As bonds saturates at doping levels above 
X ~ 0.3, which corrects the unreasonable linear decrease of the so-called pnictide height. 

PACS numbers: 74.70.Xa, 74.62.Dh, 74.62. En, 74.62.Fj, 61.05.C-, 61.43. Bn 



I. INTRODUCTION 

Superconductivity emerges in iron arsenides after the 
antiferromagnetic order of their parent compounds be- 
comes suppressed^ . This proximity of magnetic and 
superconducting order parameters is widely accepted as 
one key argument for an unconventional pairing mecha- 
nism, which is mediated by spin fluctuations^ as it was 
considered for the cuprates^^^i But in strong contrast 
to the copper oxides, superconductivity in iron arsenides 
can be induced without changing the carrier concentra- 
tion, either by applying external pressure® ^ or by charge 
neutral doping. 

Especially the isoelectronic substitution of arsenic by 
phosphorus in RE OFe{Asi^^P^) with RE = LaX^epH^ 
and ^Fe2(Asi_:rPx)2 with A = Ca, Sr, Ba, ErPHUl has 
turned out to be an ideal tuning parameter to explore 
the borderline of magnetic ordering and superconductiv- 
ity. In iron arsenides, the already present Fermi surface 
can be finely tuned by the crystal structure parameters 
in order to allocate optimal conditions for superconduc- 
tivity. However, the detailed interplay between the crys- 
tal structure, magnetic ordering and superconductivity 
is hardly understood up to now, which is to some extent 
due to the lack of precise structural data. 

One of the most intensively investigated iron-based 
material is BaFe2As2 with the tetragonal ThCr2Si2-type 
structure^. The stoichiometric compound undergoes a 
structural phase transition, associated with antiferro- 
magnetic ordering (spin-density-wave, SDW) at Ttr — 
140 K. Superconductivity emerges during the antiferro- 
magnetic ordering gets suppressed, which can be achieved 
by hole- or electron-doping- "'^®' or by applying external 
pressure^!. 



From the latter it seems easy to conceive, that also the 
so-called 'chemical pressure' by substituting arsenic for 
the smaller phosphorus atoms induces superconductivity 
in LaOFe(Asi_;j;P:r)- and BaFe2(Asi_:rPa;)^- The au- 
thors explain the appearance of superconductivity as a 
consequence of the shrinking unit cell volume in analogy 
to the effects of external pressure. On the other hand, 
the unit cell shrinks comparably by substituting barium 
for smaller strontium atoms in (Bai_xSra;yFe2As2, but 
no superconductivity appears in this casein. Thus, a 
pressure-volume effect is clearly an oversimplified expla- 
nation for superconductivity in the case of phosphorus- 
doped BaFe2(Asi_a;Pa;)2. This is the more true, as the 
crystal structures of these compounds have not yet been 
investigated in detail. We have therefore synthesized the 
series BaFe2(Asi_2,P2;)2 and (Bai_a;Sra;)Fe2As2 and de- 
termined the crystal structures by high-resolution X-ray 
powder and single crystal diffraction. We compare the 
structural data with theoretical models obtained by full- 
potential DFT calculations. We will show that the sup- 
pression of the SDW ordering as precondition to the on- 
set of superconductivity is not a simple volume effect, but 
depends on a subtle reorganization of the crystal struc- 
ture. 



II. METHODS 

Polycrystalline samples of (Bai_a;Sr2;)Fe2As2 {x = 0- 
1) were synthesized by heating stoichiometric mixtures of 
the elements (all purities > 99.9%) in alumina crucibles 
enclosed in silica tubes under an atmosphere of purified 
argon. The mixtures were slowly heated to 1123 K, kept 
at this temperature for 15 h and cooled down to room 
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FIG. 1: AC-susceptibilty of BaFe2(Asi_j;P;i,)2 with a: Ri 0.4 



temperature. The reaction products were homogenized 
in an agate mortar and annealed three to four times at 
1173 K for 15 - 25 h. The obtained black metallic pow- 
ders are slightly air sensitive and therefore handled under 
argon atmosphere. Samples of BaFe2(Asi_a;Pa;)2 with x 
= 0-1 were synthesized by solid-state reaction of the ele- 
ments likewise. The stoichiometric mixtures were slowly 
heated to 1123 K (1173 K) for 15 h. After coohng down 
to room temperature, the products were ground and an- 
nealed at 1173 K for 25 h (two to three times at 1273 K 
and 1323 K, respectively). The obtained black crys- 
talline powders show no sensitivity to air or moisture. 
EDX measurements resulted in the nominal compositions 
within 10% regarding to the Ba:Sr and As:P ratios. Sam- 
ples between x ~ 0.3 and 0.6 are superconducting with 
critical temperatures up to 29 K. No superconductivity 
was found in (B ai a; Sra; )Fe2 As2 ■ Since both phase dia- 
grams are knowii^^'^^, we show only the AC-susceptibility 
measurement of a BaFe2(Asi_:rPa;)2 sample with x sa 0.4 
as an example in Figure [l] 

The crystal structures of poly crystalline samples were 
determined by X-ray powder diffraction using a Huber 
G670 Guinier imaging plate diffractometer (Cu-Kai ra- 
diation, Ge-(lll) monochromator) . Rietveld refinements 
were performed with the TOPAS package^ using the fun- 
damental parameters approach as reflection profiles. Due 
to a better description of small peak half width and shape 
anisotropy, the approach of Le Bail and Jouanneaux^^ 
was implemented in the TOPAS program. Giving con- 
sideration to possible texture effects, the March-DoUase 
function was used for description of a preferred orienta- 
tion of the crystallites. 

Small plate-like single crystals of sa 50 x 50 x 20 mi- 
crons were selected from the polycrystalline samples and 
checked by Laue photographs using white radiation from 
a Mo-anode. Diffraction intensity data up to 29 = 
80° were collected with an Oxford Xcalibur 4-circle k,- 
diffractometer equipped with a CCD detector. Graphite- 
monochromized Mo-Ka radiation from a conventional 
sealed tube was used. The intensities were carefully cor- 
rected for absorption effects. The atom positions from 
BaFe2As2'^ were used as starting parameters and refined 
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FIG. 2: Unit cell volume against the doping levels in 
BaFe2(Asi_2,Pa;)2 and (Bai_a;Sr2,)Fe2As2. Superconduct- 
ing samples are marked by circles. Values of undoped 
BaFe2As2 under pressure are from ref!'^-' (Color online). 



by the least squares method using the Shelxl program 
packag^^*^. Positional and isotropic displacement param- 
eters of As and P were refined independently, while their 
occupation parameters were constraint to unity. 

Electronic structure calculations were performed with 
the WIEN2k package^il using density functional theory 
within the full-potential LAPW method and the general 
gradient approximation (GGA). Detailed descriptions are 
given elsewhere^Zl. LAPW is based on the muffin-tin con- 
struction with non-overlapping spheres and a plane wave 
expansion in the interstitial regions. Mixed LAPW and 
APW -I- lo (lo = local orbitals) basis sets were used to in- 
crease the efficiency of the APW linearizatiorpS. Further 
technical details are given in the monograph of SingfP^. 



III. RESULTS AND DISCUSSION 

A. SDW suppression in BaFe2(Asi_i:Pa;)2 and 
(Bai_^Sr^)Fe2As2 

The unit cell volumes of the solid solution series 
BaFe2(Asi_a;Pa;)2 and (Bai„a;Sr2:)Fe2As2 together with 
the data of the high-pressure experiments by Kimber et 
al.'23 are collected in Figure [2] Superconductivity (indi- 
cated by circles) indeed appears by pressure and by P- 
doping within the same volume range between 192 and 
199 A"^, which is marked by the dashed lines. But on the 
other hand, no superconductivity emerges by Sr-doping 
in (Bai_a;Sr2,)Fe2As2, even though the unit cell volumes 
are within the same range. Obviously, the simple shrink- 
ing of the unit cell is not the decisive condition. Su- 
perconductivity is absent in (Bai_j;Sr3;)Fe2As2, because 
Sr-doping cannot suppress the magnetic and structural 
transition. The latter becomes only shifted to higher 
temperatures, as depicted in Figure |3] 

In order to shed light on this apparent contradiction, a 
detailed look at the crystal structures is indispensable. 
First, we have refined the X-ray powder data by the 
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FIG. 3: SDW-transition temperatures of 

(Bai_2,Sri:)Fe2As2 (Color online). 
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FIG. 4: Normalized iron to pnictogen {Pn — P, As) distances 
of (Bai_a;Sri,)Fe2As2 and BaFe2(Asi_a;Pa:)2, obtained from 
Rietveld-refinements. Superconducting samples are marked 
by circles (Color online). 



Rietveld-method. The pnictogen {Pn — As, P) atoms 
occupy the 4e (OOz)-position, where zpn is the only vari- 
able parameter of the tetragonal ThCr2Si2-type structure 
(space group M/mmm). The doping dependencies of 
the normalized Fe-Pn bond lengths are plotted in Fig- 
ure |4] Sr- and P-doped BaFe2As2 show very different 
behavior, respectively. While the Fe-As bond lengths 
in (Bai_a;Sr2;)Fe2As2 remain almost constant across the 
whole doping range despite the shrinking unit cell, the 
Fe-Pn distances decrease strongly with increasing P- 
content in BaFe2(Asi_2;Pa:)2. This strongly suggests that 
a shortening of the Fe-Pn bonds is necessary to suppress 
the SDW transition. 

The latter is plausible, since it is well known that 
the Fe-As bond length drama tically influences the mag- 
netism in iron arsenides2S122l, As a first approach, we 
assume that decreasing bond lengths increase the band- 
width, thus a magnetic ground-state becomes less sta- 
ble as the bonds get shorter. Figure [5] shows sections of 
the band structures with the approximate bandwidths 
of the Fe-dxz,yz bands, which should be strongly ef- 
fected by the Fe-Pn interaction. We find that the band- 
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structures of BaFe2As2, BaFe2P2 and 
The bandwidths of the Fe-d:cz,yz bands are 



width of BaFe2P2 is much larger (+35%) compared with 
BaFe2As2. But in the case of Bao.5Sro.5Fe2As2, the band- 
width is almost identical to BaFe2As2, even though the 
volume is significantly smaller. 



Dependency of the magnetic moment on the 
Fe-As distance 



To support this assumption, we have calculated the 
dependency of the magnetic moment on the Fe-As bond 
lengths using a full-potential DFT method. It is well 
known that DFT is unable to describe the exact value 
of the magnetic moment and the correct structure at the 
same time^S. But on the other hand, it should be possible 
to get at least the changes of the magnetic moment as the 
Fe-As distance decreases. Thus we have selected a zas 
value, which gives approximately the experimental mo- 
ment of 0.8 /ie/Fe. Then we used the fixed-spin-moment 
method to calculate z-values by varying the magnetic mo- 
ments between and 1.2 /xb/Fc and relaxing the struc- 
tures subsequently. We find that the moment gets rapidly 
depleted as the Fe-As distance decreases (see Figure [6]) 
and a shrinking of 2.5 pm (« 1%) is sufficient to suppress 
the magnetism completely. This is a rather simplified ap- 
proach, since we used a ferromagnetic spin structure and 
kept the volume constant. For this reason, also the As- 
Fe-As angle varies about 1.5%, thus it is not a priori clear 
whether the bond angle or bond distance is the crucial 
parameter. 

To check the validity of our simple model, we per- 
formed full optimizations of the structure (volume, z, fi) 
using a G-type spin structure in order to model the AF 
magnetic interactions in the tetragonal phase. Figure [7] 
shows how the total energy and the magnetic moment 
depend on the bond distances and angles. The energy 
surface has a minimum somewhat beyond the experi- 
mental data (as usual for BaFe2As2), here also owing 
to the crude approximation of the magnetic structure. 
More important is the fact, that the magnetic moment 
varies strongly with the Fe-As distance, but only slightly 
with the As-Fe-As angle. This model, although still sim- 
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FIG. 6; Magnetic moments by variation of the Fe-As bond 
length in BaFe2As2 (Color online). 



to the doping value x 
Even though this 



0.25 as shown in Figure [4] 

scenario explains the difference 
between the Sr- and P-doped materials with re- 
spect to the onset of superconductivity, the data of 
BaFe2(Asi_a;P2;)2 in Figure l4] have no physical mean- 
ing. Phosphorous and arsenic cannot have the same 
coordinates in the unit cell, otherwise the atomic dis- 
tances would be unreasonable. As an example, the data 
of Figure |4] anticipate that the Fe-As and Fe-P bonds at 
X « 0.25 are about 236 pm. This is still reasonable for a 
Fe-As bond, but much too long for a Fe-P bond. On the 
other hand, both distances would be 228 pm at a; « 0.75, 
which is now acceptable for Fe-P, but impossible for an 
Fe-As bond. 
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FIG. 7: Results of the full structural optimization of 
BaFe2As2. The antiferromagnetism is approximated by a sim- 
ple G-type spin structure (Color online). 



plified, indicates that the Fe-As distance is the crucial 
parameter which determines the magnetic moment, in 
agreement with recent results reported by Johannes and 
MaziiW^. Thus, we can assume that the small shifts of 
the distance plotted in Figure [6] are associated with the 
suppression of the SDW state. 

It is worthwhile to note, that the suppression of the 
SDW in BaFe2(Asi_a;Pa;)2 may not necessarily be caused 
by the shrinking Fe-As bonds. One may also argue 
that the increasing number of Fe-P contacts enlarges the 
bandwidth, until the magnetic ground-state becomes un- 
stable. In this case, the contraction of the Fe-As bonds 
would be the consequence of the lost magnetism rather 
than its origin. However, this interpretation emphasizes 
that the Fe-As bond is in any case the gauge for the 
magnetic moment, as pointed out recently^^. 

By comparing these results with the data in Figure |4j 
we find that the shrinking of the Fe-As distances in 
(Bai_a;Sra;)Fe2As2 is only w 1 pm (-0.4%) and thus not 
sufficient to suppress the magnetism according to Fig- 
ure [6j On the other hand, the decrease of the Fe-Pn 
bond lengths is much larger in BaFe2(Asi_a;Pa;)2 and ev- 
idently sufficient to suppress the SDW. Moreover, the 
calculated value of ^ w — 1 % coincides perfectly with 
the onset of superconductivity in BaFe2(Asi_a;P2:)2 close 



Theoretical structure models 



In order to to investigate the structural behavior of 
BaFe2(Asi_a;Pa;)2 theoretically, we used crystallograph- 
ically ordered super-structures. The tetragonal space 
group PAm2 allows to calculate structural parameters for 
the doping levels x — 0.25, 0.5 and 0.75 without changing 
the symmetry. During the optimization of the atom coor- 
dinates, all compositions except a; = converged to non- 
magnetic ground-states. Thus, the SDW state is already 
suppressed at a; = 0.25 in agreement with the experi- 
ment. Furthermore, the 2;pn-coordinates of arsenic and 
phosphorus become significantly different, which lead to 
longer Fe-As and shorter Fe-P bonds, as expected from 
atom size considerations. Selected bond distances are 
collected in Table HI 



TABLE I; 


Bond lengths of BaFe; 


,(Asi_,P, 


)2 from 


DFT cal- 


culations 












X = 





0.25 


0.50 


0.75 


1 


Fe-As 


235.9 


232.7 


232.7 


232.6 


- 


Fe-P 


- 


222.0 


222.5 


221.1 


221.2 


Ba- As 


346.0 


341.2 


337.2 


335.7 


- 


Ba-P 


- 


354.5 


349.3 


342.9 


338.0 



It is evident that the Fe-As bonds shrink significantly 
only between x — Q and 0.25, accompanied by the loss of 
the magnetic moment. The contraction is about -1.4%, 
which is sufficient to suppress the SDW order as shown in 
FigurelGJ But higher doping levels do not further contract 
the Fe-As bonds, in contrast to the data in Figure |4] At 
the same time, the Fe-P bonds remain close to the val- 
ues of the pure phosphide. These findings draw a more 
realistic picture than Figure |4J because it combines the 
suppression of the SDW state by the appropriate con- 
traction of the Fe-As bond with reasonable bond lengths 
along the whole doping range. 



T3 

< 



1 


1 


•• 










1 




o 






o 


x> 





















[ J 


1 


r T 




i 






'^aTT t t t 


..... 


► A 




> 




' ^ ' J 




■< 


1 


♦ Fe-As ■■"■■■. 
■ O Fe-As calc i "~^ i^ 


o 








O Fe-Pcalc 










2 


♦ Fe-As (Sr) 











0.00 



0.25 



0.50 0.75 

x(P, Sr) 



1.00 



FIG. 8: Changes of the normahzed Fe-As and Fe-P bond 
lengths in BaFe2(Asi_j;Pa;)2 from single crystal data and in 
(Bai_a;Sra:)Fe2As2 from powder diffraction. Error bars rep- 
resent four times the estimated standard deviations. Open 
symbols represent Fe-Pn bond lengths predicted by DFT 
structure optimizations of BaFe2(Asi_a;Pi:)2. Dotted lines 
are guides to the eye. 



D. Single crystal structure determination 

From the crystal chemical considerations given in 
section B. and the results of the theoretical calcu- 
lations it is clear that the Fe-Pn bond lengths of 
BaFe2(Asi_a;Pa;)2 obtained from Rietveld refinements 
are throughout unsatisfactory, whereas the Fe-As dis- 
tances in (Bai_a;Sra;)Fe2As2 are reasonable due to only 
one arsenic position. For this reason, we have con- 
ducted high-resolution single crystal X-ray structure de- 
terminations using small crystals (tenth of microns) of 
BaFe2(Asi_^P^)2 with x sa 0.3, 0.5 and 0.7. Main 
results of the refinements together with selected bond 
lengths and angles are summarized in Table III] As ex- 
pected from chemical reasons as well as predicted by 
theory, we find the P- and As-atoms at significantly dif- 
ferent z-coordinates, resulting in different lengths of the 
Fe-As and Fe-P bonds. The initial Fe-As distance is 
240 pm and decreases rapidly up to « 25% P-doping, 
but then converges to an almost constant value around 
236 pm (-1.7%). At the same time, the Fe-P distances 
remain about 1-2 pm longer than in BaFe2P2 at least 
up to 75% P-doping. These experimental findings agree 
almost quantitatively with the theoretical prediction, al- 
though theory is not able to produce the correct abso- 
lute values. Figure IS] summarizes this reorganization of 
the crystal structures of BaFe2(Asi_a;Pa;)2 in compari- 
son with (Bai_a,Sra;)Fe2As2. A small contraction of the 
Fe-As bonds by « -1.5% (2-3 pm) is associated with the 



suppression of the SDW state, which is in turn the pre- 
condition to superconductivity. This happens by P-, but 
not by Sr-doping, where the Fe-As bonds shrink by less 
than 0.5%. 

Even though we have focused at the Fe-As bonds, we 
suggest that the Fe-P interaction is the underlying origin 
of the SDW suppression. Thus we have rather detected 
the response of the structure to the vanishing magnetic 
state. Interestingly, a similar response was not observed 
in undoped BaFe2As2 under pressure^, where supercon- 
ductivity emerges despite the Fe-As bonds remain con- 
stant as in non-superconducting (Bai_:rSr2;)Fe2As2. 



IV. CONCLUSION 

We have shown that 'chemical pressure' is an 
oversimplified view, which is unsuitable to explain 
the different behavior of charge neutrally doped 
BaFe2(Asi_a;Pa;)2 and (Bai_a;Sr^)Fe2As2. The true rea- 
son is the different response of the crystal structures to 
doping. Phosphorus doping suppresses the SDW state 
by increasing the width of the d-bands, which in turn 
leads to shorter Fe-As bonds due to its strong coupling 
to the magnetic state. No significant contraction of the 
Fe-As bond is detected in (Bai_a;Sra;)Fe2As2, where the 
SDW state persists despite the similar volume change. 
Single crystal X-ray data revealed the correct structure 
of BaFe2(Asi_2;P2:)2, where arsenic and phosphorus are 
statistically distributed, but at different coordinates zas 
and zp , which is necessary to obtain reasonable Fe-P and 
Fe-As bond lengths. The latter contract strongly at low 
doping during the magnetism is suppressed (0 < x < 
0.25) and converges to sa 236 pm at higher doping lev- 
els. This behavior is almost quantitatively reproduced 
by DFT calculations using ordered model structures of 
BaFe2(Asi_a;P2:)2. Our results strongly emphasize, that 
even subtle details of the crystal structures are crucial 
to magnetism and superconductivity in iron-based ma- 
terials, where structural data have often been obtained 
from moderate powder patterns, which allow only limited 
resolution. 
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TABLE II: Crystallographic data of 


BaFe2(Asi_a:Pi:)2 (Space 


group lA/mmm, Z = 


2, Mo-Kq radiation, A - 


= 71.073 pm) 


Compound 


BaFe2As2 


BaFe2(Aso.7iPo.29)2 


BaFe2(Aso.53Po.47)2 


BaFe2(Aso.24Po.76)2 


BaFe2P2 


Molec. mass (g/mol) 


398.88 


372.95 


357.63 


331.86 


319.77 


a (pm) 


396.24(1) 


391.78(1) 


390.65(1) 


386.60(3) 


384.35(4) 


c (pm) 


1301.35(8) 


1276.10(7) 


1273.55(6) 


1259.2(1) 


1242.2(2) 


V (nm^) 


0.2043(1) 


0.1959(1) 


0.1944(2) 


0.1892(4) 


0.1835(2) 


Abs. coefT. (mm^^) 


32.4 


29.1 


27.1 


22.5 


20.7 


Rint, Ra 


0.041, 0.024 


0.030, 0.013 


0.035, 0.016 


0.050, 0.021 


0.093, 0.042 


Ref. param., GooF 


9, 1.09 


10, 1.19 


10, 1.06 


10, 1.15 


8, 0.997 


Rl,wR2 (all data) 


0.040, 0.070 


0.023, 0.049 


0.020, 0.035 


0.033, 0.082 


0.031, 0.060 


As z position 


0.35385(7) 


0.3544(1) 


0.3542(1) 


0.3571(4) 


- 


P z position 


- 


0.3402(1) 


0.3438(4) 


0.3430(4) 


0.3459(1) 


Bond lengths (pm): 












Ba-As 


338.64(6) 


333.58(7) 


332.2(1) 


327.3(2) 


- 


Ba-P 


- 


344.0(2) 


339.6(3) 


337.4(2) 


332.4(1) 


Fe-As 


239.82(6) 


236.89(7) 


236.1(1) 


235.7(3) 


- 


Fe-P 


- 


227.2(4) 


228.9(3) 


226.0(2) 


226.14(9) 


Bond angles (deg): 












As-Fe-As 


111.40(3) 


111.57(5) 


111.6(2) 


110.2(2) 


- 


As-Fe-P 


- 


115.4(3) 


114.4(1) 


113.9(1) 


- 


P-Fe-P 


- 


119.1(3) 


117.1(2) 


117.6(1) 


116.39(7) 
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